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Abstract 

We perform a quantitative study of the neutron and electron electric dipole moments 
(EDM) in Supersymmetry, in the limit of heavy scalars. The leading contributions 
arise at two loops. We give the complete analytic result, including a new contribution 
associated with Z-Higgs exchange, which plays an important and often leading role in 
the neutron EDM. The predictions for the EDM are typically within the sensitivities 
of the next generation experiments. We also analyse the correlation between the 
electron and neutron EDM, which provides a robust test of Split Supersymmetry. 



1 Introduction 

The electric dipole moments (EDMs) of the Standard Model (SM) fermions are powerful probes of 
physics beyond the SM. Once the strong CP problem has been taken care of, the SM predictions 
for the EDMs of quarks and leptons are at least 7 orders of magnitudes below ^ the present 
experimental limits jUISlI^. The situation is drastically different in super symmetric extensions 
of the SM. The supersymmetry-breaking terms involve many new sources of CP-violation. Par- 
ticularly worrisome are the phases associated, in the universal and flavour-diagonal case, to the 
invariants ?ixg{A* Mg) and arg(^*i?). Such phases survive in the universal limit in which all the 
flavour structure originates from the SM Yukawas. If these phases are of order one, the electron 
and neutron EDMs induced at one-loop by gaugino-sfermion exchange are typically (barring ac- 
cidental cancellations [H]) a couple of orders of magnitude above the limits [01 El, a difficulty 
which is known as the supersymmetric CP problem. 

The Split limit of the MSSM |nill()|[TT] does not present a supersymmetric CP problem. Heavy 
sfermions suppress the dangerous one-loop contributions to a negligible level. Nevertheless, some 
phases survive below the sfermion mass scale and, if they do not vanish for an accidental or 



a symmetry reason, they give rise to EDMs that are safely below the experimental limits, but 
sizeable enough to be well within the sensitivity of the next generation of experiments Such 
contributions only arise at the two-loop level, since the new phases appear in the gaugino-Higgsino 
sector, which is not directly coupled to the SM fermions. 

In this paper, we perform a quantitative study of the neutron and electron EDMs in the limit 
of Split Supersymmetry. First, we compute the different contributions to the light quark and 
electron EDMs, the only relevant CP-violating operators. Indeed, quark chromoelectric dipoles 
and the gluon Weinberg operator ^2] cannot be generated at two loops. For the EDM, the 
original CP-violation in the gaugino-Higgsino sector is communicated to the SM fermions by 
gauge boson and Higgs exchanges, specifically by i) 7/1 , ii) WW, or iii) Zh exchange. No other 
possibilities are allowed at the two-loop level. 

The 7/1 exchange has been widely studied in the literature in several contexts |l,S | ll4 | n3 1 ll6j. 
The case of Split Supersymmetry was considered in ref. The WW exchange has also been 
studied in different limits jl7 | I18 | nH]. An exact 2-loop computation has been performed in the 
context of Split Supersymmetry in ref. (^Oj (see also ref. [23 ^or a computation in the context of 
a two-Higgs doublet model). Our results, for which we give explicit analytic expressions, differ 
from those in ref. [^. Moreover, we identify a third, important contribution due to Zh exchange. 
The Zh contribution is suppressed in the case of the electron EDM by a 1 — 4 sin^ 9\y factor, but 
it plays an important role in the neutron EDM. In fact, the Zh contribution is always comparable 
and often larger than the 7/1 one (which in turn is tipically larger than the WW contribution). 
We have also recomputed the QCD renormalization effect, correcting a mistake present in the 
previous literature. 



2 General expressions for the EDMs 

CP-violating phases can enter the effective Lagrangian below the sfermion mass scale rh through 
the Yukawa couplings (which are irrelevant for our study), the /^-parameter, the gaugino masses 
Mi, i = 1, 2, 3, or the Higgs-Higgsino-gaugino couplings g^, gd, g'u, g'd in 

-C = V2 (guH^W^TaHu + g'^Yn^H^BHu + gaHlW^'TaHd + g'dYn.HlBHd) + h.c, (1) 

where He = ia2H* , Ta are the SU{2) generators, and Yh^ = —Yh^^ = 1/2. The Higgs vev is in its 
usual form {H) = (0, f)^, with v ~ 174 GeV. The gaugino and Higgsino mass parameters Mi^2 
and and the couplings cju, gd-, g'm g'd are in general complex. However, only three phases are 
independent and they are associated to the invariants cpi = aig{cj'*g'^ Min) , (1)2 = arg(^*^^M2/Li), 
^ = arg{gugdgd9u)- The tree-level matching with the full theory above m gives arg(^u) = arg(^(j), 
arg{gd) = a.ig{g'^), and therefore ^ = 0, thus leaving only two independent phases. Moreover, in 
most models of supersymmetry breaking the phases of Mi and M2 are equal, in which case there 
is actually only one CP-invariant. 
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In terms of mass eigenstates, the relevant interactions are 



cw 
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where 
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(3a) 
(3b) 
(3c) 



In eq. ^^), Cf = T^f — Sy^Qj (s^ = sin^^vF) is the neutral current coupling coefficient of the 
fermion / and, accordingly, ciy± = ±cos'^9\y, Cy^+ y^- = ±(1/2 — s^)- The matrices ?7, N 

diagonalize the complex chargino and neutralino mass matrices, M_|_ = U'^M^V, Mq = N'-'^N^N, 
where Mf = Diag(M+, M^") > 0, = Diag(M°, . . . , M°) > and 



M2 (juV 

gdv /i 



/ Ml -~g'dV/^/2 ~g'^v/V2\ 

M2 gdv/V2 -guv/V2 

-g'^v/V2 gdv/V2 
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(4) 



In Split Supersymmetry, fermion EDMs are generated only at two loops, since charginos and 
neutralinos, which carry the information of CP violation, are only coupled to gauge and Higgs 
bosons. To identify all possible diagrams contributing to the EDM, let us first consider the 
case in which Mi^2,At ^ Mw- After we integrate out charginos and neutralinos at one-loop, we 
generate some effective couplings among SM bosons. These can be described in terms of gauge- 
invariant, CP-violating operators. There are 5 dimension-6 such operators: eabcW^^W^'^^Wp^ , 
H^HW^^Wt"", H'<HB^^B^"', D^H^ D^.HB'"' , Df.H^TaD^HW"'"' , where W^^ and B^^ are the 

SU{2) and U (1) gauge strengths, and and B^j^i, are their duals. The effective couplings 
relevant to generate sizable two-loop contributions to the EDM must contain 3 fields, with at 
least one photon and at most one Higgs boson. The previously-listed operators induce only the 
effective couplings 77/1, jZh, and 'yWW. Notice that CP-violating couplings of the kind 777, 
77^ and jZZ are not generated (in particular, the CP-violating operator B^^B'^^Bp identically 
vanishes unless there are three different abelian gauge fields). The absence of these couplings is 
also confirmed by an explicit one- loop calculation. Once we insert the effective couplings in a 
loop, we obtain 3 different diagrams contributing at the two-loop level to the EDM of the light 
SM fermion /, shown in Fig. ^ We therefore have 

d; = df + + (5) 



3 



^/ - Q^2^4 I™ ) TTl Iww{r^^, r^^), (6c) 

In eqs. © the sum over the indexes i,j is understood, Q/ is the charge of the fermion /, T-^fj^ is 
the third component of the weak isospin of its left-handed component. Also, rzH = {Mz / mn)'^ , 
= {M+fmnf, r+y = {M+/Mwf, r^w = {^i/^wf, where mn is the Higgs mass, and the 
loop functions are given by 

2 7o V x{l-x) J 

f I . dx . f xn + {1 - x)r2\ 

The symmetric loop function j{r,s) is defined recursively by 

rlogr j{r)-j{s) 

J[r) = -, jr,s) = . (8) 

r — 1 r — s 

In the determination of fzH we have used the symmetry of lm{DjjG^j)M^ and lm{D^jGj^)M^ 
under i ^ j. Eq. (jEJ differs from the result in ref. [20]. Analytic expressions for the functions in 
eqs. (O are given in the appendix. 

The parameters entering the expression of the quark EDM in eqs. © have to be evaluated at 
the chargino mass scale M+ . The renormalization to the scale /i at which we evaluate the neutron 
EDM matrix element is determined by the anomalous dimension of the operator qa^yj^qF^'^ . We 
find 



as{M+) 



7/26 



(9) 



where the /3-function coefficient is b = 11 — 2nq/3 and Uq is the number of effective light quarks. 
The anomalous-dimension coefficient is 7 = 8/3. To eliminate the quark mass dependence in the 
short-distance contribution, it may be preferable to consider the ratio dq/niq. Its renormalization 
is given by 

^{f^) = vi!:l^iM-'), (10) 

niq ^ TUq 

with 7' = 32/3. For as{Mz) = 0.118 it 0.004 and n = 1 GeV (the scale of the neutron mass), we 
find T/QCD = 0.75 for M+ = 1 TeV and ?7qcd = 0.77 for M+ = 200 GeV. The error in a,(Mz) 
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Figure 1: Two loop contributions to the light SM fermion EDMs. The third diagram is for a 
down- type fermion /. 



gives an uncertainty on ryqcD of about 2%, while we expect an uncertainty of about 5% from 
next-to-leading order effects. Notice that the value of t/qcd obtained here is different than what 
computed in ref. [22] and generally used in the literature. Indeed, ref. incorrectly uses the 
opposite sign for 7. Our result gives a QCD renormalization coefficient about a factor of 2 smaller 
than usually considered, and it agrees with the recent findings of ref. j23j . 

To express the neutron EDM in terms of the quark EDMs, we use the results of QCD sum-rule 
techniques [2il l25]: 



where /,r ~ 92 MeV and we have neglected the contribution of the quark chromoelectric dipoles, 
which does not arise at the two- loop level in the heavy-squark mass limit. Note that dn depends 
on the light quark masses only through the ratio mu/md, for which we take the value niu/md = 
0.553 ± 0.043. 



3 Expansions in the heavy-chargino limit 

We now discuss the result in the limit in which the i?-symmetry breaking scale, determining 
gaugino and Higgsino masses, is larger than Mz and niH- A leading-order perturbative expansion 
ofeq. © in powers of |Mi_2/^|/-^l and |Mi^2/u|/m^ (keeping all orders in |Mi^2//^| and in Mz/ruff) 




(11) 
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gives 

f2 



" ^^^^4^ il^'^'^^'^^^l^^'T^'T^^rj 



where gu,d-, Qud-' -^1,2 and /U now indicate the absolute value of the corresponding quantity and 
the functions F^h-, Fzh, Fww are given in the appendix. As long as M1/M2 < 1 (as, for 
instance, in the case of gaugino masses unifying at the GUT scale), the second term in eq. (|12b ) 
is suppressed with respect to the first term and numerically it is not very significant. Notice that 
eqs. (|12|) explicitly exhibit the dependence on the two CP-violating invariants \gugd/M2fJ.\ sin (p2 
and \g'^gy Ml fj,\ sin (pi. Because of the suppression of the second term in eq. (fT^ l. both the 
electron and neutron EDM are mostly characterized by a single invariant. 

While eqs. (jT^.b) can be obtained from an expansion at the first order in v/Mi^2, eq. H12b ) 
arises only at the second order. This is because the origin of df^ can be traced back to the 
vertices W^^ Dfj^HDyH ^ B^'^D^HD^H in the unbroken electroweak symmetry phase, requiring 
two insertions of the Higgs vev. The additional M/v factor in eq. (|Ht) leads to a contribution 
to the EDM which is parametrically of the same order of dj^ and d^^ in the v/M expansion. 
Notice also that the coefficients of the potentially large logarithms of |Mi^2At|/"i/f correspond 
to the anomalous dimensions that mix the EDM operator with the CP-violating dimension-6 
operators obtained from integrating out the supersymmetric particles well above the weak scale. 

The relative importance of the three contributions to df in eq. © can be estimated from 
eqs. (|12|). Let us consider for definitess the case M2 = /i. By keeping only terms enhanced by a 

large log(M2/i/m^) in the expressions (PT|) for F^h, Fzh, F^\y, we obtain Fzh ~ F-yni—Syy + 
(2) 

3/4) and F^-^^ « —F^h/^- As a consequence, we find 



^ (T3,,-2.^Q,)(3-4.^^) dj^ ^ _ 

where in the expression for dj^ we have neglected the subleading second term in eq. (|12b '). 

Numerically, eq. ^ gives d^^ « 0.05 d2^, df^ ^ -0.3 dl^ and df^ ^ dZ^ , d^^ « -0.7 dZ". 
These simple estimates show the importance of the ZH contribution to the neutron EDM. A 
detailed numerical analysis of the relative importance of the different contributions to the elec- 
tron and neutron EDMs is given in the next Section. The qualitative estimates above are in a 
remarkably good agreement in the large M2 = /x limit. 

4 Numerical analysis 

We now perform a numerical analysis of the full results for the EDM in eqs. We consider 
a standard unified framework for the gaugino masses at the GUT scale. By using the RGEs 
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given in refs. |1U1 126j . the parameters in eqs. © can be expressed in terms of the single phase 
(p = (p2 and the four following positive parameters: M2, // (evaluated at the low-energy scale), 
tan/5, and the sfermion mass scale rh. In first approximation, the dipoles depend on (3 and 
(f> through an overall factor sin 2/J sin Therefore, in order to maximize the effect, we choose 
tan/3 = 1 and the phase (p evaluated at the low-energy scale such that sincf) = 1. Notice that 
if the ratio |/i/M2| is much larger or much smaller than one, the renormalization effects which 
mixes /x and M2 (peculiar of Split Supersymmetry) tend to suppress the effective phase. In this 
case, a maximal CP violation can only be achieved with very particular choices of the initial 
values for the higgsino and gaugino masses. Once we have fixed sin 2/3 sin (/> = 1, we are then 
left with the three dimensionful parameters M2, fi, rh. The overall sfermion scale rh enters only 
logarithmically through the RGE equations for gu^d, d'ud- choose to present the results as 
contour plots in the M2-^ plane and set rh = 10^ GeV, which is consistent with the cosmological 
bounds given in ref. [27j . 

Fig. 121 shows the prediction for the electron EDM, the neutron EDM, and their ratio dn/de- 
The red thick line corresponds to the present experimental limits d^ < 1.6 x 10~^^ecm while 
the limit dn < 0.63 x 10~^^ecm |2] does not pose a constraint on the parameters shown in Fig. |2j 
In the Split limit, and assuming gaugino mass unification, all EDMs are controlled by a single 
phase. The results for de and dn shown in Fig. |2l scale approximately linearly with sin 2/3 sin 0. 

A robust test of Split Supersymmetry can be performed if both the electron and the neutron 
EDM are measured. Indeed, in the ratio dn/de the dependence on sin(/>, tan/3 and m approx- 
imately cancels out. This can be easily understood from eqs. (jl2l) which show that, as long 
as the chargino and neutralino masses are sufficiently larger than Mz, the only dependence of 
dn/de on gu4, g'^^, and sin 0^2 comes from the existence of the Mi-dependent term in 
see eq. p2b ). This term, as previously discussed, is numerically small. Nevertheless, because of 
the different loop functions associated to the different contributions, the ratio dn/de varies by a 
O (100%) factor when the M2 and fi are varied in the range spanned in the Figures. Still, the 
variation of dn/de is comparable with the theoretical uncertainty on the determination of dn in 
terms of quark EDMs in eq. (|1H) due to the hadronic matrix element, and is significantly smaller 
than the variation in the ordinary MSSM prediction, even in the case of universal phases (see 
e.g. ref. |2H]^)- On the other hand, the usual tight correlation between the electron and muon 
EDMs, dfj^/de = m^/rrie persists. 

Fig.EJshows the relative importance of the different contributions to the EDMs. For the same 
reason explained above, the results shown in Fig.|31are fairly independent of sin (p, tan /3 and rh. As 
anticipated, the ZH contribution to the electron EDM is suppressed by the T-^f^ — 2s'^Qf factor. 
On the other hand, the corresponding contribution to dn is always important and represents the 
largest contribution in a significant portion of the parameter space shown in Fig. El The WW 
contribution is also sizable, especially in the case of the neutron EDM, but is typically smaller 
than the ZH or ones. While the ZH and ^H contributions always add constructively, the 
WW contribution has an opposite sign. However, also due to the ZH contribution, its size is not 
large enough to flip the sign of the overall electron or neutron EDM. 



^Note that the ZH contribution is missing in the analysis of the Spht Supersymmetry case in ref. |28|. which 
leads to a stronger correlation between de and dn. 
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Figure 2: Prediction for d„, d^, and their ratio dn/dg. We have chosen tan/3 = 1, sin(/> = 1, and 
m = 10^ GeV. The results for dn and scale approximately linearly with sin 2/3 sin while the 
ratio is fairly independent of tan/?, sin(^ and rh. The red thick line corresponds to the present 
experimental limit < 1.6 x 10~^^ecm jl]. 
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Figure 3: Relative importance of the different contributions to the EDMs. We have chosen 
tan/3 = 1, siiKp = 1 and m = 10^ GeV, but the result depends only very weakly on this choice. 
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5 Conclusions 



In summary, we performed a quantitative study of the neutron and electron EDMs in Split Su- 
persymmetry. Clearly, our results also apply to the MSSM in any limit in which the contributions 
involving sfermion exchange are suppressed. Our result for the WW exchange differs from the 
one in the literature. Moreover, we find a new contribution associated to Zh exchange, which 
plays an important and often leading role in the neutron EDM. We have also given the correct 
value for the QCD renormalization of the quark EDM. We performed an analytical and numerical 
analysis of our results, summarized in Figs. [21 and 01 The correlation between the electron and 
neutron EDMs is found to be stronger than in standard super symmetric scenarios, and it may 
become a crucial experimental test for Split Supersymmetry. Still, we find an O (100%) variation 
of dn/de in the parameter space we considered, which is comparable with the hadronic uncertainty 
on the determination of dn in terms of quark EDMs. The results summarized in Fig. [2 are quite 
promising in the light of the expected impressive improvement of the experimental sensitivities 
in the years to come and represent one of the most relevant windows on Split Supersymmetry. 
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Appendix 

We write below the analytic expression of the loop functions in eqs. © and Q. 
1 



fiH{r) 



VI -4r 



log r log . ■ , - + Li2 



1 - VI - 4r 



Li2 



1 + Vl-4r, 



(14) 



fzH{r, ri,r2) 



r - 1 
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logr 
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log 1 



Xi 



log 1 
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(15) 



fwwiri,r2) 
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log r2 



log 1 
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log 1 
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1 — Xl 
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where 



r — ri + r2 =F Y^(r — ri + r2)^ — 4rr2 
yi,2 = r2 + (n - r2)xi,2 



— 

Vi = yi\r=l- 



(17) 

(18a) 
(18b) 



In the physical region, the functions /y//, fzH and fww do not develop imaginary parts. 

In order to study the case of chargino masses larger than Mz , or ri , r2 ^ 1 , it is useful to 
switch to the variables R = ^/rir^, P = '>'i/r2 and expand in 1/R. We then get 



2R 



i?2 



fzH{r,ri,r2) = azH[p)^- + ^ + 

logi? 



/M^VK(n,?^2) = aww{p)—;:r- H 7^ V O 



R 



R 



R^ 



where 



aznip) 
bzH{r, p) 
aww{p) 
bww{p) 



^^P 



2(p-l) 
(p-l)2 
(p-l)2 



logp 

r log r log /9 
1 -r 

(p - 1 - logp) 
(p+1 



/O-I + 



Li2(l-p)+Li2(l-l/p) 



log p + Li2(l - p) - Li2(l - 1/p) 



(19a) 
(19b) 
(19c) 

(20a) 
(20b) 
(20c) 
(20d) 



F^u{p.R) 

FzH{r, p, R) = Azh{p) log R + BzH{r, p) + O 
FwwiP^ R) = ^wwip) log R + B^D^ip) + O 



FivUp^ R) = ^I^V(P) log R + ^w(p) + O 



logi? 

R 
\ogR 

R 
logR 
R 



From these expansions, the expressions for the functions F^jj, Fzh-, PwW' ^ww ™ P^P 
follow: 

logii' 



,(2) 



(2) 



(21a) 
(21b) 
(21c) 
(21d) 
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where 
Azh{p) 
BzH{r, p) 



(P- 1)(2 - p) - plogp 



+ 



4(p-l)2 

^ {(2-2r + rlogr)(p-l) [p 2 - 2s2^,(p _ 1)] 



(1) 



A. 



4(r-l)(p 
+ (r-l)(p-l) [| + l-52^(p+l) 

+ (r-l)p[Li2(l-p)-Li2(l-l/p)]} 
-p2 + l + 2plogp 



log p + rp log r log p 



(2) 



8(p-l)3 

-4p(p - 1) + (p2 - 4p - 1) logp - 4p [Li2(l - p) - Li2(l - 1/p)] 

16(p- 1)3 
(p-7)(p-l) + 2(p + 2)logp 



wvr/-r 8(p-l)3 

4(p - 4) (p - 1) - (p2 + 4p + 7) logp - 4(p + 2) [Li2(l - p) - Li2(l - 1/p)] 



P 



16(p- 1)3 
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(22b) 

(22c) 
(22d) 
(22e) 
(22f) 
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